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INTRODUCTION 


Under Task II o£ contract NAS 6*^2173 WOLF was respon** 
sible for assisting NASA/Wallops Flight Center in test plan- 
ning, data reduction and analysis for the Marine Navigation 
Systems Evaluation Experiment (MANSEE3 . Under a .previous con- 
tract, NAS 6-1942, WOLF developed the mathematical techniques 
to reduce and analyze the MANSEE data. See Appendix 1 for a 
mathematical description of this program. The results of this 
analysis are summarized in this report. 

Section 1.0 describes the MANSEE test. Section 2.0 
describes the navigation sensors which were exercised. 

Section 3.0 describes the ground truth instrumentation and 
the processing of ground truth data. Section 4.0 summarizes 
the residual statistics for individual navigation sensors. 

Here residuals were calculated by differencing the actual 
measurements with anticipated measurements computed from the 
ground truth trajectory. In most cases the ground truth 
trajectory was sufficiently accurate that these residuals 
could be attributed to navigation sensor measurement errors. 

Section 5.0 presents some of the results obtained by 
using the SEAMAP program to filter data from the navigation 
sensors. The resultant filtered trajectories were differenced 
with the corresponding ground truth trajectories to obtain 
navigation position and velocity errors. Section 6.0 con- 
tains various conclusions and recommendations derived from 
the MANSEE effort. 

Volume 2 contains the results from SINS data. This 
was made a separate volume because the results are classified. 
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SECTION 1.0 

TEST DESCRIPTION AND CONDENSATION 


The purpose o£ the MANSEE experiment was to exercise 
and evaluate representative navigation systems in current 
general usage by the maritime community. The test vehicle 
was the Apollo Range Ship USNS VANGUARD, Available naviga- 
tion systems and navigation aid subsystems on-board the 
VANGUARD were supplemented with rented or borrowed naviga- 
tion instrumentation so that representative systems from 
each genre were exercised. The test range was NASA's 
Wallops Flight Center, located at Wallops Island, Virginia, 
where precision tracking radars are available on the coast 
to provide the necessary "Ground Truth," or nominal trajectory 
of the ship (see Figure lA) . 

The test was performed in three phases during June 
1972. Phases one and three were instrument calibrations 
conducted at Port Canaveral, Florida, the home port of the 
VANGUARD,, Phase two was the experiment operation, conducted 
on June 5, 6, and 7, 1972, in the immediate vicinity of 
Wallops Island, Virginia. 

The design goal of the experiment procedure was to 
simultaneously exercise the various navigation systems on 
the ship through a series of varying trajectory geometries 
and velocities. To maintain ground- truth control, the ship 
remained within line of sight of the land-based radars, a 
distance of 20-25 miles. To avoid shallow water, the ship 
remained 10-15 miles from shore. Thus, a sector of ocean 
approximately 15 miles by 25 miles was available for opera- 
tions. The ship executed rectangular, crossing, and circular 
patterns at various velocities within the sector (see Figure 
IB). 
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The design goal of obtaining simultaneous data from 
all available sensors was met during several long periods 
of time during the test, with the exception of the relative 
velocity sensor. The EM Log was inoperable during the experi- 
ment. The available data set at the conclusion of the test 
contained in excess of one million independent position fixes 
for the three day operation - an amount more than adequate 
for sensor evaluation. 

The performance of the navigation sensors was within 
or close to the nominal accuracy for all systems except 
Loran-C and the AN/SRN-9, Apparently the Z- channel of 
Loran-C locked up on the wrong cycle of the carrier pro- 
ducing a 2500 meter error. However, the Y-channel was 
accurate to within 24 meters (RSS) . The usable data from 
the AN/SRN-9 consisted of only four points which is insuf- 
ficient to state definite conclusions of its accuracy. 


f 


5 



SECTION 2.0 
INSTRUMENTATION 



2.1 NAVIGATION SENSORS 

The selection o£ the VANGUARD as a test vehicle was 
advantageous both from the standpoint of having a highly 
trained crew to operate the navigation systems, and having 
a complement of navigation sensors available, as the follow- 
ing list indicates: 

1. LORAN-C Receiver 

2. Inertial Navigation System (SINS) 

3. Marine Star Tracker 

4. Electro-Magnetic Velocity Log 

5. Mark- 19 Gyrocompass 

6. SRN-9 Navy Navigation Satellite System 
Receiver 

To this instrumentation set was added, through lease 
contracts or subcontracts, or loans from other government 
agencies, the following navigation sensors: 

1. Raydist Navigator 

2. Motorola Mini -Ranger 

3. Omega Receiver 

Detailed specifications on »:he above instrumentation 
are available from Wallops Station, A summary description 
of each system follows: 

LORAN-C - The LORAN-C Receiver on the VANGUARD is a 
Collins SPN-3S. 

t 
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SINS 


The inertial navigation system on the 
VANGUARD is a Sperry MK3 Mod S. The 
system consists of four gyros and three 
accelerometers mounted on a stabilized 
platform, all of which are input to a 
MINDAC computer. The output from the 
MINDAC software consists of both gyro 
torqueing signals and the ship^s position 
(latitude and longitude) , velocities 
(north, east, vertical, and total hori- 
zontal) and heading, 

MST - The Marine Star Tracker is a subsystem 
of the SINS. The instrument measures 
azimuth and elevation to stars for input 
to the MINDAC computer. The MINDAC soft- 
ware outputs ship's position (latitude 
and longitude) based on the star measure- 
ments . 

EM LOG - This is a sensor on the VANGUARD for mea- 
suring ship's velocity relative to the 
surrounding water. The EM LOG was not 
operating during the MANSEE experiment. 

MARK-19 GYROCOMPASS - The MARK-19 is a Sperry instru- 
ment. Measurements are ship's heading 
relative to the earth's spin-axis (true 
north), plus ship’s roll and pitch angles 
relative to the local gravity vector. 

SRN-9 RECEIVER - The AN/SRN-9 receiver on the VANGUARD 
was manufactured by ITT. The system receives 
and decodes transmissions from the Navy 
Navigation Satellite System. The measurements 
and satellite ephemeris are subsequently 
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processed (together with velocity and 
heading from the INS or the MK 19/EM LOG) 
in the UNIVAC 1230 computer on-board the 
ship. Output from the computer is a 
position fix (latitude and longitude) and 
standard deviation of the fix. 

RAYDIST - The RAYDIST >'T" navigator supplied by 

TELEDYNE HASTINGS - RAYDIST was used in 
the experiment. The navigator is a phase 
measuring receiver and was used in conjunc 
tion with the lower Chesapeake network of 
RAYDIST base stations. Measurements are 
transformed to latitude and longitude by 
determining the intersections of lines 
of position. 

MINI -RANGER - The Motorola Mini-Ranger is a line-o£- 
sight two-way range measuring system. The 
configuration used in this experiment had 
the range console, receiver/ transmitter , 
and antenna on-board the VANGUARD, with 
three transponders mounted at surveyed 
locations on land. 

OMEGA - The OMEGA receiver used in .the experiment 
was a Dynel Electronics Model 200 loaned 
to NASA by the U , S . Navy , 


2.2 "GROUND TRUTH" INSTRUMENTATION 

"Ground Truth" instrumentation consisted of two 
AN/FPS-16 C-Band pulsed tracking radars, and two C-Band 
transponders. These radars, normally used to track rockets 
during launch and satellites in orbit, have been qualified 


8 



as geodetic instruments, with measurement accuracies o£ 

2-3 meters in range and iS-30 arcseconds in angles. One 
radar is a permanent element of the VANGUARD'S instrumenta- 
tion complex, and a C-Band transponder was mounted on a 
tower on Wallops Island as a point source target for the 
radar to track. The other system is part of the V/allops 
Island radar complex, and a C-Band transponder was mounted 
on the vanguard's superstructure as a target for this radar. 
Detailed specifications on this instrumentation and its 
utilization are available from Wallops Flight Center. 


> 
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SECTION 3.0 

GROUND TRUTH DATA PROCESSING AND ANALYSIS 


The data was subsampled by merging in time sequence 
alternate one hour segments of data on computer tapes. An 
analysis of all the ground truth data at a one per second 
rate on these tapes was conducted to validate the consistency 
of the ground truth measurements. The ^ priori error statis- 
tics (measurement noise and bias) for the radars from the 
pre-mission and post-mission calibrations were compared to 
the mission data. 

The shipboard azimuth measurements contained a simis^oi- 
dal error having peak variations of 20 to 30 mils as compared 
to systematic errors of less than .1 mil for the land based 
radar. A probable explanation is that the ship's heading 
must be determined to a high degree of accuracy to process 
the measurement. Because of the wide discrepancies in data 
quality it was decided to use for ground truth data processing 
only shipboard range measurements and land based azimuth 
measurements. The error statistics used in all ground truth 
processing are as follows: 

Bias 

0 .5 meters 

0 .1 rails 

The propagation of these errors along with possible survey 
uncertainties of +_ 6 meters in latitude and longitude of 
the position of the land transponder yields a ground truth 
accuracy of the ship's trajectory of 5.2 meters (one sigma) 
relative to the position of the land radar. The ship's 
velocity' is determined to + .033 meters/second (one sigma) 

— t 

and the ship's heading to + 1.0 degree (one sigma). 


Noise 

RADAR RANGE (Ship) + 2 meters RMS 

RADAR AZIMUTH (Land) + .2 mils RMS 
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Figure 3.1 shows the ship trajectories with respect 
to the Wallops FPS-16 and Mini-Rangers. Figure 3.2 shows 
the RMS Cpi’edictedj position error vs. time for the 1422 
data span when processing the radar data. The transponder 
location error one sigma contribution is 3.5 meters and the 
one sigma contribution due to measurement errors (noise and 
bias) is about 3.9 meters. Figure 3,3 shows the correspond- 
ing RMS velocity error vs, time. 

Table 3.1 lists the radar residual statistics for 
all the data spans. Note that Wallops FPS-16 range (iv'hich 
was not weighted in the solution) agreed closely with the 
range measured by the Vangard FPS-16. The only exception 
was the 1612-1645 data span which had a mean range residual 
of 7.8 meters. It is not known why the bias was so large 
in this run since the bias was fairly constant throughout 
the run. 
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■ Figure' .3,2. 
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Table 3.1, Radar Residual Statistics For Individual Data Spans 

(Meters and Mils) 


Data Span 

No. Data 
Points 

Wallops FPS 

Range Range 
Mean Sigma 

-16 

Azimuth 

Mean 

Azimuth 

Sigma 

No, Data 
Points 

Vangard FPS' 

Range Range 
Mean Sigma 

-16 

Azimuth 

Mean 

Azimuth 

Sigma 

1422-1500Z 

219 

-0.1 

2.0 

-.01 

0.25 

211 

0.0 

1.9 

1.93 

12,74 

1612-1645Z 

152 

7.8 

3.0 

0,00 

0.45 

1 

145 

0.09 

2.0 

22.03 

24.30 

1912-1940Z 

ISO 

0.4 

1.8 

-0.03 

0.29 ■ 

141 

-0.1 

1.9 

-4.48 

16.43 : 

2100-2200Z 

351 

0.1 

2.3 

-0.01 

0.57 

■351 

j 

0.0 

1.9 

-3.66 

15.33 

1133-1212Z 
(June 7) 

225 

3.4 

4.1 

-0.03 

0.74 

225 

0.6 

3.1 

37.82 

44.74 

OVERALL 

1097 

1.8 

2.8 

-0.01 

0.52 

1073 

0.1 

2.2 ■ 

9.50 

25.35 


SECTION 4.0 

NAVIGATION SENSOR RESIDUALS 


IVhen processing multiple sensor data with the SEAMAP 
program, measurements from any sensor may be either used by 
the filter to improve the state estimate or processed in a 
residual analysis mode. For a ground truth run, data from 
both C-Band radars is used to update the state estimate while 
data from the navigation sensors is treated in the residual 
analysis mode. 


A detailed residual analysis of the navigation sensors 


was performed for the 

following 

seven data 

June 

6 

0030Z 

- 0130Z 

June 

6 

1422Z 

- 1500Z 

June 

6 

1612Z 

- 1645Z 

June 

6 

1S46Z 

- 1900Z 

June 

6 

1912Z 

- 2000Z 

June 

6 

2100Z 

- 2200Z 

June 

7 

1133Z 

- 1212Z 


During these periods good quality ground truth data was avail- 
able and most of the navigation sensors (except the EM Log) 
were operating. 


• While performing the data reduction., it was discovered 
that the location of the C-Band, radar transponder on board the 
ship was not input to the SEAMAP program. Since the trans- 
ponder was 50.6 meters from the reference point of the ship, 
the computed ship position was in error by about 50 meters. 
Unfortunately, the data from the 0030-0130Z and 1846-1900Z 
spans were lost before the runs could be remade. The effect 
of this error on the results was computed (approximately) from 
the ship's heading, but results from these data spans are not 
listed here except for sensors which have errors much larger 
than 50 meters. The antenna locations of all other sensors 
were included in SEAMAP in all of the runs. 
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4.1 


MEASUREMENT RESIDUAL MODEL 


For each Navigation sensor, the measurement residual 
is computed as 



1.^ CXj .13 


( 4 . 1 . 1 ) 


where : 

Actual measurement vector 

State vector estimated from C-Band radar data 

a priori values of sensor systematic error 
parameters . 

If we express the measurement residual in terms of its com- 
ponent errors and neglect higher order terms: 




A, 6s ■- M.6X. 
1 ” J “1 


( 4 . 1 . 2 ) 


whe7 5 : 


r^* = Random measurement error having zero mean 

and covariance matrix 

6£ = Error in a priori values of sensor systematic 

error parameters having zero mean and 
covariance matrix S 



Error in state vector having zero mean and 
covariance matrix 
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= Partial derivatives o£ the measurement with 
respect to the systematic error parameters 
in s 


Mj = Partial derivatives of the measurement with 
respect to the adjusted parameters in Xj • 


IVhen a sensor is treated in the residual analysis mode, these 
three error vectors are statistically independent. 


In running SEAMAP , it is the usual practice to initiate 
residual analysis after the filter has converged and is 
nearly constant. 

A 

For the residuals analyzed, and 6Xj are neandy 

constant over a typical data span. In modeling the sensor 
errors, it is assumed that ^ is constant. Thus, each residual 
is hypothesized to be the sum of a random error plus two 
error terms that are constant over the data span. However, 
these constant terms may vary between data spans. 
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4.2 


RESIDUAL STATISTICS 


To determine the magnitude of the error terms, for 
each data span, the mean and standard deviation of each 
component of the residual vector were computed as 



(4.2.1) 


where n^|^ is the number of residuals in the i 


. th 


data span 


(i=l,2,..7) and 5y ■ - represents a component of the residual 


ij 


vector at time t^ during the i 


, th 


data span. 


The overall mean and standard deviation for each 
residual component were computed as 



(4.2.2) 


where : 


n 


total number of data points 
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i=l 
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Note that the overall standard deviation contains two com- 
ponents; one due to the scatter of residuals about the mean 
within the data span, and the second component due to the 
variability of the data span means about the overall mean. 
This second term is due to variations in the systematic 
residual components between data spans. 

From equation 4.1.2 the covariance matrix associated 
with the mean residual vector for a given data span is 
given by 


ECu.u. 3 = V. = MPM^ + ASA^ + — R C4.2.3) 

^—1—1 1 n . ^ 

For the data spans analyzed, n^ was sufficiently 
large that the measurement noise contribution — could be 

m X 

neglected. The state error contribution MPM"^ is analytically 
computed within SEAMAP. The presence of significant systematic 
measurement errors can be assessed by computing the ratio 


g 


u. 

1 




n,n 


C4.2.4) 


where u. is a residual mean and q is the corresponding 
X ri j n 

T 

diagonal element of the matrix MPM . 

For large n- , if the moan measurement error is zero, 
the statistic "g" follows a normal distribution with mean 
zero and standard deviation one. Thus, if the magnitude of 
the residual mean is much larger than /q^ ^ , v/e can say 
that the sensor bias is significantp.y different from zero 
and the value of \x^ is the best available estimate of the 
corresponding element of A6^. The standard deviation of 
this estimate is approximately /q^^^ , 
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Under the assumptions that the error vectors 6X. and 
A ”J 

s and the arrays Aj and Mj are nearly constant over a data 
span, the expected value of the residual standard deviation 
is the corresponding measurement noise standard deviation. 
These assumptions are clearly violated ivhen the residual 
plots show significant trending or when the residual noise 
characteristics change during a data span. 


I 
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Table 4.1 Navigation Sensors Overall 
Measurement Residual Statistics 


Sensor 

Measurement 

Mean 

Sigma 

S amp 1 e 
Size 

MINIRANGER 1 

Range 

"2.8 meters 

3.6 meters 

62 

MINIRANGER 2 

Range 

-1,5 meters 

4.4 meters 

217 

MINIRANGER 3 

Range 

-5.0 meters 

5 . 3 meters 

552 


Channel Y 

-21 meters 

12 meters 

605 

LORAN-C 

Channel Z 

-2581 meters 

133 meters 

956 

MK-19 GYRO 

Heading 

3,3 degrees 

2,3 degrees 

1421 


Latitude 

100 meters 

551 meters 

4 

AN/SRN-9 

Longitude 

-170 meters 

575 meters 

4 


Stations B-C 

-1040 meters 

611 meters 

107 

OMEGA 

Stations D-F 

262 meters 

589 meters 

107 


Red Network 

23.8 meters* 

2.9 meters 

1632 

RAYDIST -T 

Green Network 

-11.1 meters* 

3.5 meters 

1631 


*For Raydist, the mean residuals are primarily due to 
initialization error rather than system errors. 


t 
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FIGURE 4.1. TYPICAL SEAMAP RESIDUAL PLOTS DAY 15B - 1422 DATA SPAN 





4.3 


MINI-RAWGER RESIDUALS 


IVhen considering the results listed in Table 4.1, 
Table 4,2 and Figure 4.2, particular emphasis should be 
placed on the ranges designated Mini“Ranger #1. The 
transponder for this range link v/as mounted on the same 
tower as the C-Eand transponder for the ship radar, pro- 
viding essentially colinear range measureme-its . Thus the 
data could be compared without possible errors due to 
survey, geometry, or other measurements. The result yields 
a high degree of confidence in this range link. There is 
a definite possibility that the errors observed in the 
other range links (see also Figure 4.2) are due to a source 
such as survey errors, and that their error characteristics 
are actually similar to #1. If both Mini ^2 and Mini ^3 
had approximately the same bias, there would have been 
reason to suspect an azimuth bias on the IVallops FPS-16. 
Since the biases were of opposite signs, an azimuth bias 
on the land radar seems unlikely. It is also possible that 
the errors observed were due to transponder antenna height 
or pointing. Motorola cautioned that transponder antenna 
pointing could be critical and excessive height could cause 
problems due to multipath. Furthermore, the Mini- Rangers 
were operating at the limit of their range. 
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Table 4.2, Mini-ranger Residual Statistics For Individual Data Spans 


Data Span 

Mini #1 

No. Data 
Points 

(Meters) 

Mean Sigma : 

Mini # 2 

No. Data 
Points 

(Meters) 

Mean Sigma 

Mini #3 

No. Data 
Points 

(Meters) 

Mean Sigma 

0030-0129Z 

— 

— 

1 

— 

— 

— 

— 

— 

— 

1422-lSOOZ 

— 

— 


167 

-5.8 

3.4 

82 

-0.5 

8.1 

1612-1645Z 

— 

— 

— 

3S 

11.8 

6.2 

194 

-5.2 

4.6 

1912-2000Z 

S3 

-3,8 

3,6 

7 

18.5 

4.5 

218 

-6.2 

4.9 

2100-2200Z 

— 

— 

— 

— 

— 

— 

— 

— 

— 

1133-1212Z 

9 

3,0 

3.6 

8 

13.8 

5.0 

58 

-6.3 

4.4 

(June 7) 










OVERALL- 

62 

-2.8 

3.6 

217 

-1.5 

4 , 4 

552 

-5.0 

5.3 






ASUREMENT RESIDUAL (METERS} 



4.4 


LORAN-C RESIDUALS 


Figure 4.3 and Table 4.3 show the mean and one sigma 
spread o£ the LORAN-C residuals over all data spans. The 
0030 span was eliminated from the Channel Y display 
(Figure 4.3) because of the large transients in this 
channel. From this plot, it is evident that the systematic 
errors in both channels did not change significantly between' , 
data spans. 

The values listed in Table 4.3 include the secondary 
phase corrections and ASF corrections (see Appendix, Section 
3,5). Therefore, the large Channel Z residual (-2581 m) is 
probably due to the receiver locking up on the wrong cycle. 
Since the carrier frequency is 100 KHz, one cycle is 10 ysec. 
In the MANSEE test area, 10 ysec delay from Channel Z rep- 
resents 2431 meters. Had the receiver been properly locked, 
the true error would have been about 150 meters. Even this 
is larger than should be expected and may be due to an 
additional calibration error in the receiver since the 
results are quite consistent. It is also possible that the 
ASF correction (which accounts for propagation over land 
instead of seawater) was slightly in error. This was 
supplied by the U.S. Coast Guard and is computed as a 
"best fit" to measured data. For Channel Z, the slave 
station is located in Dana, Indiana, and the master station 
is i^ iarolina Beach, Worth Carolina. Thus, both the signal 
from the master to the slave and from the slave to the ship 
must travel over 1000 km of irregular terrain (Allegheny 
Mountains) which has varying conductivity. Furthermore, 
the path from the master to the ship is along the coast 
which also has varying conductivity. However, in spite 
of these problems, the error in the ASF correction should 
be less than 0,2 ysec or 49 meters. 

The standard deviation of th^ Z-channel residuals is an 
order of magnitude larger than the Y residual sigma. This is 
a result of the frequent -200 m transients in the Z-channel 
which may also be related to improper receiver calibration, 
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Note the consistency of Channel Y (slave is on 
Nantucket Island) . The mean residual only varied from 
“17 to -30 meters with 6 to 17 meters noise. For this 
channel, the scale factor in the MANSUE test area is 
ISO m/ysec. Thus, the -21 meter overall mean is an error 
of only 0. 14 ysec. 
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Table 4.3, Loran-C Residual Statistics For Individual Data Spans 



Channel 

Y (Meters) I 

Channel 

2 (Meters) 


No. Data 
Points 

Mean 

Sigma 

No. Data 
Points 

Mean 

Sigma 

OOSO-OISOZ 

— 

— 

— 

348 

-2570 

181 

1422-15002 

179 

-18 

10 

179 

-2616 

148 

1612-16452 

59* 

-20 

17 

60 

-2579 

52 

' 1846-19002 

80 

-30 

9 

80 

-2565 

58 

1912-2000Z 

! 

96* 

-17 

6 

97 

- 2580 

58 

2100-2200Z 

191* 

-21 

14 

192 

-2575 

71 

1133-12122 
fJune 7) 

— 

— 

— 

— 

— 

— 

OVERALL 

60S 

-21 

12 

956 

-2581 

133 


* 


One bad Channel Y point was deleted 




FIGURE 4.3. LORAN-C MEASUREMENT RESIDUALS 
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MK-19 RESIDUALS 


From Tabic 4.1 the overall mean MK-19 heading residual 
is 3.3 degrees and the standard deviation is 2,3 degrees. 

From Table 4.4, the mean residuals for each data span varied 
from 2.1 degrees to 6 degrees and the corresponding standard 
deviations varied from 1.8 degrees to 3,6 degrees. When 
compared with the ground truth heading accuracy of 1 
degree Cone sigma), it can be concluded that the MK-19 
heading measurements contained a positive bias. Since 
the MK-19 measures ship heading, while the ground truth 
computes the velocity vector, any ocean current would 
result in a bias between the two. However, the bias was 
fairly constant in all the data spans even though the 
ship changed headings. The residual plots showed no 
evidence of trending (Figure 4..4) . Therefore, we can 
conclude that the residua!! standard deviations within 
each data span are good estimates of the measurement 
noise sigma. 


f 
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Table 4.4. MK-19 Gyrocompass Heading Residual 

Statistics For Individual Data Spans 



HEADING DEGREES 

DATA 



NO. DATA 

SPAN 

MEAN 

SIGMA 

POINTS 

0030“0130Z 

2.7 

2,2 

348 

1422-1500Z 

6.0 

1.9 

211 

1612-1645Z 

2.4 

3. 6 

145 

1912-2000Z 

4.7 

2.0 

141 

2100-2200Z 

4.6 

1.8 

351 

1133-1212Z 
(June 7) 

2,1 

2.8 

225 

OVERALL 

3.3 

2.3 

1421 
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4.6 AN/SRN-9 RESIDUALS 


The AN/SRN-9 provided psoition fixes at approximately 
one hour intervals during the experiment. Unfortunately, 
June 6th was the only day on which simultaneous ground truth 
and AN/SRN-9 fixes were available. The mean and standard 
deviation of four position residuals are given in Table 4.1 
and plots of residuals vs. time are given in Figure 4.5. 

The mean latitude residual is 100 meters and the standard 
deviation is 551 meters. The mean longitude residual of 
-170 meters can be attributed partly to the different geoids 
used within the AN/SRN-9 (NWL-8D) and within the SHAMAP 
Program (FISHER), This longitude difference on the North 
American continent is 51.4 meters at the MANSEE exercise 
latitude Thus, after adjusting for geoid differences, 
the actual mean longitude residual is 128.6 meters. Due 
to the limited number of data points, it is difficult to 
determine if the AN/SRN-9 errors are primarily random or 
systematic . 


TTl 

■^Anderle and Smith, *'NWL-8 Geodetic Parameters Based On 
Doppler Satellite Observations,” NIVL Report 2106, July 
1967. 


‘-tJ' 
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FIGURE 4.5. 
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OMEGA RESIDUALS 


Table 4.5 lists the OMEGA residual statistics for all 
data spans and Figure 4.6 is a plot of the residuals. It 
is seen that the random error is generally smaller than the 
systematic component. The largest source of error is due 
to errors in the propagation correction. From the plots 
of the propagation correction versus time (see Appendix) , 
it is seen that the correction changes drastically during 
the transition from night to day over the propagation path. 
Thus, the large errors during the 0030-0130Z and 1133-1212Z 
spans are due to errors in the rapidly changing propagation 
correction. However, it is interesting to note that the 
worst case errors are still within the stated accuracy of 
1 nautical mile during the day and 2 nautical miles at 
night, despite the fact that the ship is only 350 nautical 
miles from station D (Forrestport , New York),* 


The propagation corrections are somewhat less accurate for 
ranges less than 450 nautical miles. 



Table 4.5 


. Omega Residual Statistics For Individual Data Spans 
(Including Propagation Corrections) 


Data Span 

Stations B- 

No. Data 
Points Mean 

-C 

Sigma 

Stations D- 

No. Data 
Points Mean 

F 

Sigjua 

0030-0130Z* 

23 

"117 

746 

23 

3067 

946 

1422-1500Z 

16 

-557 

997 

16 

-696 

99 

1612-1645Z 

10 

736 

339 

10 

-309 

176 

1S46-1900Z* 

6 

-1287 

109 

6 

-1085 

120 

1912-2000Z 

11 

-1457 

156 

11 

-1298 

114 

2100-2200Z 

25 

-1459 

308 

25 

-295 

209 

1133-1212Z 

16 

-2926 

207 

16 

-9 

344 

OVERALL 

107 

-1040 

611 

107 

262 

« 

589 


^Ground truth is in error by about 50 meters. 
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4.8 


RAYDIST-T RESIDUALS 


Table 4.6 and Figure 4.7 show the mean and one sigma 
spread of the Raydist residuals. Considering that the MANSEE 
test area was outside the prime coverage area of the Raydist 
network (see Figure 4.8), the random Raydist residuals are 
quite small (Approximately 0,03 lanes). Since the ground 
truth, RMS position error due to measurement noise (radars) 
and state noise is about 3 meters, it cannot be stated whether 
the Raydist residuals are due to errors in the ground truth or 
to Raydist error. The nominal accuracy of the Raydist network 
is 0.02 lanes which is about 2 meters in the MANSEE test area.* 

The spread of the mean residuals was 8*7 to 34,2 meters 
for the Red Channel and -1.8 to -16.7 meters for the Green 
Channel. The spread of the means is far more meaningful than 
the actual values of the means since the Raydist receiver was 
not accurately initialized during the test. Thus the overall 
means of 23.8 meters (Red Channel) and -11.1 (Green Channel) 
are simply indications of the initialization error. The spread 
of the means is somev/hat larger than expected. In particular, 
the 8.7 meter Red Channel mean for the 1612 data span and the 
“16,7 Green Channel mean for the 1422 data span stand out. 

Both of these mean residuals are difficult to explain. For 
the 1612 data span, the Red Channel measurement is indicating 
that the ship should be about 16 meters northwest of the ground 
truth position. This is the only data span in which the 
Wallops FPS-16 range did not agree with the Vanguard FPS-16 
range. However, if the Wallops radar is correct, the actual 
ship position is 8 meters southeast of the ground truth position. 


*In the MANSEE test area, one lane is 112.0 and 82,1 meters 
respectively for the Red and Green Channels, 
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Likewise, in the 1422 data span, the Raydist Green 
Channel indicates that the ship is about 10 meters north of 
the ground truth position, Miniranger "2 tends to confirm 
this, but the accuracy of the Minirangers is not sufficient 
to state a definite conclusion, 

I 

During the MANSEE exercise, the Raydist receiver 
temporarily lost track on two occasions and had to be reset 
manually by examining strip chart recordings. The first time 
was on the evening of June 5 and the second was early morning 
on June 7. Thus the manual resetting could not explain the 
spread of the mean Raydist residuals. There was also a period 
after 1400Z on June 6 when the Days Point station (Green Channel) 
was transmitting at low power. This might explain the discrepancy 
in the Green Channel mean for the 1422 data span. Apparently 
the receiver lost track between 1422Z and 1612Z because both 
channels have a +1 lane offset between the 1422 data spans 
and the other data spans. Likewise, there was a -1 lane, 

Green Channel offset in the 2100Z data span. 


( 
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Table 4*6. Raydist Residual Statistics For Individual Data Spins 



Red 

Channel 


Green Channel 

Data Span 

No. Data 
Points 

Mean 

Sigma 

No. Data 
Points 

Mean 

Sigma 

1422-1500Z 

729 

24, 2 

3.0 

728 

"16.7 

2.5 

1612-1645Z 

253 

8.7 

3.0 

253 

-7.8 

5.6 

■1912-200CZ 

235 

28.9 

1.7 

235 

-1.8 

3.1 

2100-2200Z 

200 

24.0 

1.5 

200 

-6.7 

2.0 

1135-1158Z 

215 

34.2 

4.1 

215 

-10.3 

4.3 

OVERALL 

1632 

23.8 

2.9 

1631 

-11.1 

3.5 
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SECTION 5.0 

NAVIGATION FILTER RESULTS 


The previous section gives detailed measurement 
error characteristics for individual sensors. This section 
summarizes the overall navigation accuracy achievable with 
various sensors and combinations of sensors. Here the 
sequential filter in SEAl'lAP was used to process a time 
sequence of navigation measurements and to give a current 
estimate of ship position and velocity Cx>y!)» Error 

analysis computations within the program converted sensor 
measurement error statistics into predicted navigation 
error levels (RJIS position error, RMS velocity error). 

The actual navigation errors position error, RSS 

velocity error) were evaluated by differencing the estimated 
trajectory with the corresponding ground truth trajectory. 


5.1 COMPARISON OF ACTUAL AND PREDICTED NAVIGATION ERRORS 

Figures 5.1, 5.2 and 5.3 show comparisons of actual 
position errors (RSS) and predicted error levels (RMS) , 

In filtering the Mini-Ranger data (Figures 5.1 and 5.2), 
the RIIS error was of the same magnitude as the RSS error 
when two Mini -Rangers were tracking and their error char- 
acteristics were correctly modeled. The Ri^IS error was 
computed assuming 4 meters random noise, 3 meters bias, 

4 meters latitude error and 3.2 meters longitude error. 

In Figure 5.1, the RSS error was more than double the RMS 
error during the period from 1442 to 14492. During this 
period, the tracking from Mini-Ranger “3 became intermittent 
and its error characteristics suddenly changed (see Figure 
4.1). The bias appeared to change by almost 10 meters after 
1442Z. Likewise in Figure 5.2, the jRSS error suddenly grew 
when Mini- Ranger j/ 1 dropped out and Mini-Ranger It 2 began 
tracking. Mini-Ranger 112 had a bias of +19 meters for this 
span (Table 4.2). 
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In filtering the Loran-C data, all propagation cor- 
rections were applied to the data and the one cycle CIO 
microseconds) error in the Z channel was corrected so that 
the results more realistically model actual horan performance. 
Figure 5.3 shows the comparison of RSS and RMS position errors 
assuming 1 ysec random noise and 0.2 ysec phase correction 
error for each channel (other error sources were . negligible) . 
The RSS errors ranged from 104 to 207 meters and during most 
of the span, they were approximately 50% greater than the 
predicted (RMS) error. The major source of error is the 
150 meter, uncorrected bias in the Z channel. As discussed 
in Section 4.4, this represents about 0.6 ysec error which 
is three times larger than the propagation error assumed in 
computing the RMS position error. Hence, the discrepancy 
between R^^S and RSS errors is expected. 


I 
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POSITION ERROR (METERS) 
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FIGURE 5.1. MINI-RANGER FILTER RMS AND RSS POSITION ERRORS - 1422 DATA SPAN 
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FIGURE 5.2. MINI-RANGER FILTER RMS AND RSS POSITION ERRORS - 1912 DATA SPAN 
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5.2 COMBINATION OF VELOCITY DATA WITH LORAN-C POSITION 

DATA 

Figures 5.4, 5.5 and 5.6 show the effect of using the 
SEA^^AP filter to combine the Loran position data with velocity 
data accurate to 0.1 m/sec. These results represent an upper 
bound on the navigation accuracy improvement derived from 
combining velocity data with Loraii“C position data. 

The 1618 data span (Figure 5.4) is somewhat unusual in 
that the RSS error when using velocity data is greater than 
with Loran data only. The second Loran-C, Z channel, data 
point in this span was in error by 10 ysec (2431 m) . Since 
this occurred at the very beginning of the run when the data 
editor was turned off,* the position estimate was grossly 
in error. As additional Loran data was processed, the 
position estimate returned close to the correct position. 
However, the inclusion of velocity data made the filter 
overconfident in its estimate and the filter gain matrix 
became far too small. Thus new Loran position fixes were 
essentially ignored. 

Figures 5.5 and 5.6 are more typical of the improve- 
ment expected when velocity data is added to Loran position 
fixes. The improvement is not great because the dominant 
error is the ISO meter bias in the Z channel of Loran. In 
fact, Figure 5.6 shows that the RSS error, when using 
velocity data, can exceed the Loran only error for a con- 
siderable period of time. This is again due to the reduction 
in the filter gain matrix when the additional measurement 
type is added. 


The data editor is inoperative for the first 100 seconds of 
the run so that initialization errqrs do not result in the 
rejection of good data. 
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FIGURE 5.4. LORAN-CAND LORAN-C + VELOCITY DATA 
FILTER POSITION ERRORS, 1618 DATA SPAN 
lALL CORRECTIONS APPLIED TO DATA) 
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FIGURE 5.5. LORAN-CAND LORAN-C + VELOCITY DATA 
FILTER POSITION ERRORS, 1912 DATA SPAN 
(ALL CORRECTIONS APPLIED TO DATA) 
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FIGURE 5.6, LORAN-CAND LORAN-C + VELOCITY DATA 
FILTER POSITION ERRORS, 2100 DATA SPAN 
(ALL CORRECTIONS APPLIED TO DATA) 
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DATA EDITING EXAMPLE 


SEAMAP contains a statistical data editor which 
compares each measurement residual vector with its analyti 
cally computed covariance matrix. Outlying components of 
the residual are deleted prior to processing by the sequen 
tial filter. During the 1612 data span an equipment mal- 
function in the LORAN receiver introduced a -100 ysec 
(-15,000 m) bias in the Channel Y data (see Figure 5.7). 
This data transient caused the filter position estimate 
to be in error by 20,000 meters, With the data editor 
activated, all of the bad Channel Y data was deleted and 
the filter solution, using only Channel Z data, showed 
no noticeable deterioration during the twenty minutes 
following the initial transient. The inclusion of 
velocity data would have allowed the filter to follow 
reasonable changes in ship velocity during this period 
of bad LORAN data. 
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FIGURE 5.7. PROCESSING OF LORAN--C 1G12TIME SPAN SHOWING EFFECT OF DATA EDITOR 5S 


SECTION 6.0 
CONCLUSIONS 


In performing the MANSEE experiment and subsequent 
data analysis, NASA/Wallops station has demonstrated the 
test planning, instrumentation, and data processing 
capability to utilize multiple sensors in establishing a 
ground truth ship trajectory accurate to _+ 5 meters (one 
sigma). This value is considered to be a "worst case" 
one sigma, as the survey uncertainty of six meters used 
in the error propagation is a pessimistic value. The 
accurac^ was achieved by using a generalized sequential 
filter (SEA^^AP) to process combined data from the two 
C-Band tracking radars. The availability of higher data 
rates with the Wallops radars indicates the potential 
for obtaining the same or better accuracy for numerous 
other navigation and traffic control testing applications, 
particularly in aeronautics. 

Following the experiment, data analysis with the 
SEAMAP program demonstrated the desirability of redundant 
ground truth instrumentation. Comparison of azimuth 
residuals w'ith pre and post-mission calibration data 
showed that the ship radar azimuth measurements contained 
systematic errors two orders of magnitude greater than 
those for the land based radar. Elimination of ship azimuth 
measurements reduced ground truth position errors from 10-20 
meters IIMS to 5 meters RMS. 

The usefulness of a sequential filter in processing 
measurements subject to intermittent data losses or transients 
was demonstrated. In the case of LORAN-<", SEAMAP successfully 
processed through 20 minutes of bad Channel Y data \vhile 
providing an adequate navigation solution. 
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The scope o£ this effort did not include an indepth 
analysis of all data from all sensors, but rather attempted 
to indicate the overall characteristics of each. Thus, the 
results presented here only summarize the operational error 
characteristics of the particular navigation instrumentation 
installed on the VANGUARD during the MANSEE test. 

Most of the navigation sensors performed as expected. 
The only exceptions were the Z channel of Loran-C (which 
apparently locked up on the wrong cycle of the carrier) 
and the AN/SRN-9 (which had only 4 data points). To sum- 
marize, a comparison of the overall RSS residuals obtained 
from ground truth processing with the nominal accuracy of 
the systems is given below: 


Nominal Accuracy 


RSS Residual 


Mini-Rangers 

Loran-C 

MK-19 GYRO 

AN/SRi>l-9 

Omega 

Raydist-T 


3 m 

15-400 m 

7 

100-370 m 
1852 m 
2 m 


6,5 m 
2023 m 
4.0“ 
580 m 
967 m 
3,2 m* 


In the case of the Mini-Rangers and Raydist-T, the 
RSS residual cannot be interpreted as the true error because 
of 5 meter uncertainty in the ground truth. 


Does not include mean error which is'due to improper initializa 
tion. 
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SECTION 1.0 
INTRODUCTION 


The ’’SEAMAP" program is designed to process multi- 
sensor maritime navigation data to obtain an optimal esti- 
mate o£ ship position and velocity 'i'be program 

has extensive error analysis capability for optimal or sub- 
optimail treatment of sensor static error parameters such as 
measurement bias and station location errors. 

This report contains a mathematical description of 
the sequential filter, sensor models and state dynamic model 
contained in "SEAMAP." 

Figure 1.1 shows overall structure of a sequential 

filtering program which processes an ordered time sequence 

of measurements J ^ ••• 'J’be state dynamic 

model computes the state and the transition matrix at time 

t given the state at t . . From this state the sensor 
n n- 1 

models compute the a priori or anticipated measurement and 
the partial derivatives of this measurement with respect to 
the state. The filter combines these quantities v/ith the 
actual measurement to obtain an updated or aposteriori esti- 
mate of system state. 

In "SEAMAP" the filter, dynamic model and each sensor 
model have been programmed as independent software modules 
such that any one can be modified or replaced without dis- 
turbing the remaining code. 
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SECTION 2.0 

SEQUENTIAL FILTER DESCRIPTION 


If navigation filter accuracy is predicted under the 
assumption that sensor errors are entirely due to zero mean 
random noise previous WOLF studies (2) have indicated that 
the presence of additional static errors, such as sensor 
bias, can result in considerable deviation of the actual 
filter accuracy from the predicted values. In many cases 
these static parameters are not observable and therefore 
cannot be included in the state vector of adjusted parame- 
ters. Johnson (Ij has derived a generalized linear filter 
which considers a priori knowledge of these unobservable 
parameters in computing an optimal sequence of filter weight 
matrices . 

These previous developments have been incorporated 
by WOLF into the design of SEAMAP . 

SEAMAP contains a generalized linear filter which 
processes multisensor marine navigation data. The filter 
has extensive capability for treating static error para- 
meters in the sensor models (such as measurement bias and 
station location errors) . Each static error parameter may 
be treated by the filter as: 

• Adjusted - the static parameter is included in 
the state vector of parameters to be adjusted. 

• Consider - a considerable parameter is not 
adjusted but its uncertainty is taken into 
account where computing the filter weighting 
matrices. This means that the filter is able 
to make optimal use of any a priori knowledge 
about the uncertainty in parameters which are 
unobservable . 
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© Analyze Only - error analysis computations 

external to the filter compute the additional 
uncertainty in each adjusted parameter due to 
the uncertainty in each analyze only parameter. 

The display of these uncertainties (one sigma) 
gives the detailed numerical consequences of 
omitting a parameter from the filtering compu- 
tations. 

Classification of each static parameter is designated 
by card input for any given filter run. 

Outlying data points may cause the filter to diverge 
or go out of control, SEAMAP contains statistical data 
editing logic which deletes outlying measurements. 

The follov . ■ sections describe the generalized linear 
filter equations . implemented in SEAMAP. In these discus- 
sions and flow charts the subscripts of each of the vectors 
or matrices refer to first the data point at which the matrix 
or vector is being evaluated, and second, the data point at 
which observation data was last incorporated into the estimate. 
Thus means the estimate prediction of X at data point i, 

based on observation data through data point i-1. 
means the estimate of X at data point i-1, based on observ'’a- 
tion data through data point i-1. §. . is evaluated at 

1 5 1 “ i 

X. 1 = X. 1 T . 

1-1 i-l/i-l 


i 
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2.1 


STATE AND MEASUREMENT MODELS 


The state vector can be partitioned as: 



where : 

= dynamic force model terms such as position and 
velocity which must always be included in the 
set of adjusted parameters. 

= adjusted and consider parameters such as 
measurement bias and station location errors 
for sensor number one. 

2 

= same for sensor number two. 


Xj^ = same for sensor number i. 


The elements of X^ must always be adjusted while the 
adjustment of any static error parameter is optional. 

It is assumed that the incremented state model has 
the form: 


X*^ = f (X'^ „ O + q 
— n — n n n-l-^ 


where g_ is a zero mean vector of state noise having covariance 

matrix Q . 

n 
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A measurement from the k 
current state by: 


th 


sensor is related to the 







+ 



k 


where r. is a zero mean noise vector having covariance 
matrix and ^ contains those static error parameters 
of the kth sensor which are treated as unadjusted analyze. 


2.2 FILTER EQUATIONS 

For a typical filter cycle it is assumed that a priori 
errors in the adjusted parameters can be approximated by a 
linear combination of four statistically independent error 
vectors. 


6X 


■n/n-1 ^n/n-1— ^ ^'^n/n-l^^/o ^ “U/n-l ^ — n/n-1 


where : 


^n/n-1— ~ component due to errors in unmodeled 

analyze static error parameters. 

W n 6X / == component due to errors in the initial 

n/n-i — 0/0 ^ 

state . 

= component due to zero mean state noise. 

^n/n-1 ~ component due to zero mean measurement 

noise. . 
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Under this assumption, the total a priori state covariance 
matrix is given by 


^^n/n “ ^n/n-1 ^ ^n/n-1 ^'^n/n-l ^o/o ^'^n/n-1 


^ ^n/n-1 ^ *^n/n-l 


In SEAMAP the last 3 components of grouped 


as 


^n/n-1 “ ^'^n/n-1 ^o/o ^'^n/n-l ^ ^n/n-1 ^ ^n/n-1 

At time t„ the force model subroutines fs^e flow chart 
Figure (2.1)) evaluate 2^.n/n-l ^n n-1 previously 

estimated state 5Ln_i/;n-l* quantities are then used to 

update the unmodeled analyze coefficient matrix. 

^n/n-1 %,n-l ^n-l/n-1 

The state covariance components due to initialization 
errors, state noise and sensor noise are updated as 

P / i=® iP 1/ 

n/n-1 n,n-l n-l/n-1 n,n-l ^n 


The sensor model subroutines then evaluate the measure- 
ment residual, 5y.j^> partial derivatives of the measurement 

with respect to the adjusted parameters, M , the partial 
derivatives of the measurement with respect to unadjusted 
analyze parameters, A^, and the measurement noise covariance 
matrix R„ . 



The residual covariance matrix is computed as 



E 


C6y 6y^) = M P . 

^ — n-^ n n/n 


M 


n 


R 

n 


A chi square statistic, c^, is computed to test the 
consistency of the residual, 6y , with its analytically com- 
puted covariance matrix, C , If e is above some preset 
^ n n 

threshold, the editing logic is activated; otherwise, the 
measurement residual is used to adjust the state vector. 

The adjust gain is computed as: 


T - 1 

K = P 

n n/n-1 n n 


The rows of corresponding to consider parameters 
are filled with zeroes and the adjustment is given by 


—n/n —n/n - 1 




Using the above gain matrix (with consider rows zeroed) the 
unadjusted analyze coefficients are updated as: 


= (I-K 11 ) L . T - K A 
n/n ^ n n n/n-l n n 

And the after adjustment state covariance matrix is computed 
as 


Vn ' Pn/n-lfl-W 


K R K ' 
n n n 
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L 


! 


f 


For display purposes the total state covariance matrix, 
including the effects of uncertainties in the unadjusted 
analyze parameters is computed as; 


V / = p, I , slT 

n/n n/n n/n n/n-1 
2.3 DATA EDITING LOGIC 

Each measurement residual vector is compared with the 
analytically computed residual covariance matrix by computing 
the statistic 


c 


n 




Under the assumptions used to derive the sequential filter 
follows a chi square distribution with degrees of freedom 
equal to the number of elements in the measurement vector 

This statistic is scale free (i.e., independent of the 
measurement units such as meters, seconds, degrees, etc.) and 
the significance levels depend only on the number of degrees 
of freedom. 

A very large value of may be due to 

1. A single outlying measurement. 

2. Use of incorrect noise statistics within the 
measurement model. 

An unmodoled state distuprbance . 
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3. 



% t t i I I h . i 

The data rate was very high Cl/second) in comparison 
to the response time of the Vanguard to unmodeled forces. 
Because of the redundancy of ground truth quality sensors 
it was possible to verify the measurement error statistics 
for each sensor. Therefore in reducing the MANSEE data it 
was assumed that large values of were only caused by 
outlying measurements. 

Whenever a large value of is detected the compo- 
nents of 6};^^ are compared with their standard deviations. 
Those components outside the ^3 sigma points are deleted 
from the current measurement. The threshold value is input 
for each sensor. 


( 
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State vector of adjusted parameters 



s 





— i 



Vector of sensor model static error parameters 
Cassumed zero mean) 

Measurement vector 

Anticipated measurement vector 

Nonlinear vector function giving state vector at 
time i as a function of state vector at some 
earlier time 

Nonlinear vector function giving the anticipated 
measurement as a function of the current state, 
and of the sensor static error parameter ^ 

State transition matrix computed in the force 
model subroutines 


= Sensitivities of the adjusted parameters wrt the 
sensor static error parameters (s) 


S = ECs"^ s) 



Total state covariance matrix 
State noise covariance matrix 


Sensor noise covariance matrix 



Measurement* residual 




Partial derivatives of the measurement wrt the 
vector of adjusted parameters 


A. = Partial derivatives of the measurement wrt the 

1 

unadjusted sensor static error parameters 


£ fjv 

Qo 


Residual covariance 
Filter gain matrix 
Residual chi square 
Rejection threshold 
Nominal state noise 


matrix 

statistic 

for current sensor 
covariance matrix 


Component of state error due to zero mean 
noise 




state 



Component of state covariance matrix due to zero 
mean state noise 


G 


i/j 


Component of state covariance matrix due to zero 
mean state noise 


P 


i/j 


State covariance matrix not including effects 
of unmodeled analyze parameters. 


* 
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SECTION 3.0 
MEASUREMENT MODELS 


SEAMAP currently includes measurement models for the 


following 

sensor types: 

1 . 

Ship velocity (x,y) 

2. 

Ship position 

3. 

Land based radar 

4. 

Ship based radar (R,A^) 

5. 

Range only (R) 

6. 

LORAN-C (At^, At 2 ) 

7. 

Hyperbolic Lane Count AJI 2 ) 

8. 

Ship heading 


The program is designed so that additional sensor 
types can be included vvithout disturbing any of the existing 
code . 


As described in Section 2.0 each sensor measurement 
is related to the system state by the nonlinear vector 
equation; 




1 ) 




where 

= force model parameters \\'hich are ship latitude, 
longitude, east velocity and north velocity. 

= static parameters of the sensor which are 

—1 ^ 

to be adjusted or treated as consider paia- 
meters . 
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Ic t ll 

s - static error parameters of the k sensor 

which are treated as unmodeled analyze 

parameters . 


k 

r. = zero mean measurement noise vector. 

The sensor subroutines compute a priori values of and 

k d ^ 

the partial derivatives of with respect to 

and s^. 

— 1 — 


3.1 SHIP VELOCITY MEASUREMENT MODEL 


SINS measures ship east and north velocity. In 
addition to random errors it is assumed that each measure- 
ment contains an offset error plus a bias that varies linearly 
with time. Thus the a priori velocity measurements are 
modeled as: 


^1 

= X + + 


At 

" ’^1 

^2 

- y + 


At 

* 


where : 




Ui,a2 




r 


1 



At 


actual east and north velocities, 
east and north offset errors, 
linear bias coefficient, 
zero mean random noise terms, 
time since filter ini*tialization . 
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3.2 


SHIP POSITION MEASUREMENT MODEL 


Direct measurements of ship latitude and longitude 
are obtained from SINS, AN/SRN-9 and the Marine STARTRACKER. 

In addition to random noise it is assumed that each measure- 
ment contains an offset error plus a bias that varies linearly 
with time. The a priori measurements are modeled as: 


>^1 = 

^ + otji^ ’*■ 


At 

* ^1 


X + tt2 + 

h 

At 

* ^2 


V 

where : 


= actual latitude and longitude. 

defined as above. 

3.3 LAND BASED RADAR MODEL 

The land based radar measures range and azimuth to a 
ship mounted transponder. Each of these measurements contains 
a constant bias plus random noise. In addition the station 
location is subject to survey errors. Thus the radar measure- 
ments are given by: 


+ AR + r^ 

+ AA^ + 1*2 
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where . 


<j),X = ship latitude, longitude. 

= station latitude, longitude. 

^ A(}> ,AX = station location survey errors. 

s s 

ARjAA^ = Range, azimuth bias errors. 

^2 ” random measurement noise. 

Range, azimuth and the appropriate partial derivatives are 
computed within subroutine ARC. 


3.4 SHIP BASED RADAR 

The ship based radar measures range and azimuth rela- 
tive to the ships bow o£ a transponder on shore. These 
measurements are assumed to contain a constant bias plus 
noise. The location o£ the transponder is subject to survey 
errors. Thus the a priori ship radar measurement is computed 
as: 


= Rf(}),X,({)^,X^,A(}i.p,AX^) + AR + r^ 

Y2 = AC(j),X,(i),p,X^, A(t)^,AX^) - ATAN (x/y) + AA^ + 



where : 


= transponder latitude, longitude 

“ transponder location survey error 

X - ship east velocity 

y - ship north velocity 

and the remaining symbols are defined as above. 

Range, azimuth and the appropriate partial derivatives 
are computed within subroutine ARC. 


3.5 LORAN-C 


LORAIvI-C is a pulsed hyperbolic positioning system. 

The measurements are the differences in arrival time between 
a master pulse and the pulses from each of two slave stations. 
The station locations are subject to survey errors. Each 
time difference measurement contains a bias terra which is 
the sum of secondary phase correction errors and interstation 
delay errors. The index of refraction is also treated as a 
systematic error parameter. 


The LORAN-C measurements are modeled as: 


n 


“ c ^^s^ ~ ^m ^ ^ ^"1 


“ E„ + Ef m D + r. 


s, -m "s^m “1 


^2 = I C's - \ * E - E_ . E, „ . D_. . r. 


'm s,m ^ms. 


) 
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where 


c - velocity of light in vacuum 

n - index of refraction (n=l. 000338) 

R ,R = ranges from ship to slave //I, slave #2 
^1 ^2 

R = range from ship to master station 

m D JT 

R ,R^ ^ranges from master to slave “1, slave #2 
in s ^ jti s 2 

E ,E secondary phase corrections for slave "1, slave i>2 

Sf s^ 

D ,D =fixed delay between reception of master pu-^se 

IT13 <1 m S 

and transmission of slave pulse for slave nl, 
slave ^f2. 

r^,T 2 = zero mean measurement noise. 

The values for E ,E_ ,E .E ,E are computed analytically 
assuming propagation over seawater'- . 

R<100 statute miles (161 Km) 

E = 821.544/R - .011402 + 1.0936x10"^ R 

R^lOO statute miles (161 Km) 

E = 38673. 4/R - 0.40758 + 2.15476x10"^ R 

where R is in meters. 


^^^Winkler, G ,M. R. /’Path Delay, Its Variations, and Some 

Implications for the Field Use of Precise Frequency Standards," 
Proceedings of the IEEE, Vol, 60,»No. 5, May 1972. 


A- 18 



In addition to the secondary phase correction E, 
there is another correction which takes into account the 
fact that propagation may be over land rather than over 
seawater. This correction (called the AST) must be computed 
for each pair of stations at th ^ particular position of 
interest. Values obtained from the U.S. Coast Guard for 
the MANSEE test area are: 

Channel Y: ASF = -1.229 ysec 

Channel Z; ASF = +2,088 ysec 

For simplicity, these numbers are combined with the fixed 
delay between stations so that the values used for D are: 

Channel Y: D = 49000.-1.229=48998.771 ys 

smi 

Channel Z: D = 65000 , +2 . 088=65002 . 088 ys 

sm^ 

These measurements and the appropriate partial derivatives 
are computed within subroutine LORP . 

3.6 HYPERBOLIC LANE COUNT MEASUREMENTS 

Hyperbolic lane count measurements are obtained from 
OMEGA and RAYDIST. These systems are similar to LORAN-C 
except that the basic measurement is the phase difference 
between two CW signals rather than the arrival time difference 
for two pulses, A lane count measurement is ambiguous in that 
the measuring receiver only establishes position within a lane 
(One lane = A/2 meters where X is the continuous wave length.) 
Thus the actual measurement is some fraction of a lane. The 
integer lane count is arbitrarily set at the beginning of a 
navigation exercise and is updated either manually or auto- 
matically each time and the receiver records a lane crossing. 
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SEAMAP obtains the initial integer lane count from 
an input parameter card. As a default option this integer 
count is taken to be the integer part of the first lane 
count measurement processed. OMEGA and RAYDIST station 
locations are subject to survey error. Lane count measure- 
ments contain bias errors which are the sum of the second- 
ary corre'~tion errors for each station. The phase velocity 
is also treated as a systematic error parameter. 

Within SEAMAP, a priori lane count measurements are 
computed by first computing the arrival time differences 
and then converting these differences to lane counts. Thus 
the lane count is given by: 




rz = ^ - ''Bj^ " \ - \ * ^2 * ^2 


ivhere : 


V 




£ 





r 


1»^'2 





phase velocity 

distances from ship to channel 1 stations 

distances from ship to channel 2 stations 

CW transmission frequency 

Channel 1 secondary phase corrections 

Channel 2 secondary phase corrections 

zero mean measuremefit noise 

lane count origin adjustment 
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For OMEGA, v was assumed to be 300574 000 and £ \v'as 10,3 KIlz. 

The predicted propagation corrections (E^ ,E^ ,Ejj ,Eg ) for 

12 12 

OMEGA are tabulated by the Defense Mapping Agency Hydrographic 
Center, Figures A-1 and A-2 are plots o£ the propagation 
correction for the MANSEE test area during early June 1972 
for stations B,C,D and F. At the time of the test, station D 
was located in Forrest Port, New York and station F was in 
North Dakota. Both stations were operating on temporary 
assignments . 

The secondary phase corrections for Raydist were 
assumed to be zero due to the shorter distances involved. 

V was assumed to be 299670000 and £ was 3.307400 MHz (Red 
Channel^ and 3.307550 Mliz (Green Channel). These a priori 

•■U 

measurements and the appropriate partial derivatives are 
computed within subroutine MELCNT. 


3.7 HEADING MEASUREMENTS 

Ship heading relative to north is measured by SINS and 
the MK- 19 gyrocompass. In addition to random errors it is 
assumed that heading measurements contain an offset plus a 
bias that varies linearly with time. Thus the a priori 
heading measurement is given as : 

= ATAN (x/y) + a + B At -*• r 

where : 

x,y = actual east, north velocities 
a, 3 = offset and time linear term 

X = random measurement noise 

This measurement and the appropriate partial derivatives are 
computed witliin subroutine MEIIEAD. 
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